A Fully Quantum Mechanical Model of a SQUID Ring Coupled 
to an Electromagnetic Field 
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A quantum system comprising of a monochromatic electromagnetic field coupled to a SQUID ring 
with sinusoidal non-linearity, is studied. A magnetostatic flux <& x is also threading the SQUID ring, 
and is used to control the coupling between the two systems. It is shown that for special values of 
the system is strongly coupled. The time evolution of the system is studied. It is shown that 
exchange of energy takes place between the two modes and that the system becomes entangled. 
A second quasi-classical model that treats the electromagnetic field classically is also studied. A 
comparison between the fully quantum mechanical model with the electromagnetic field initially in 
a coherent state and the quasi-classical model, is made. 



I. INTRODUCTION 

With the SQUID ring (here taken to be a thick super- 
conducting ring enclosing a single Josephson weak link 
device) regarded as having potential for future quantum 
technologiesii2i2ii, it is clearly of interest to consider its 
interaction with an external quantum mechanical electro- 
magnetic (em) field. This interest has certainly been pro- 
moted by the recent experimental work on the creation 
of quantum mechanical superposition states of Josephson 
systems£*§iL&, with particular emphasis on the existence 
of such states in SQUID ring o 9 i 10 . As with these latter 
experiments, in order to investigate these states we con- 
sider a monochromatic em field with frequency u> e (typ- 
ically in the 0.1 to 1 THz region), coupled to a SQUID 
ring oscillator with frequency u> s . In addition a magne- 
tostatic flux <& x is also applied to the ring, as depicted in 
figure 1. Since the primary purpose of the work reported 
here is to study the full quantum mechanics of this cou- 
pled system, we make the assumption that the operating 
temperature (T) is such that hw e 3> fcgT, huj s ^> ksT 
so that both the ring and field modes behave quantum 
mechanically. 

As we shall show, in this fully quantum mechanical de- 
scription the quantum states of the em field mode plus 
SQUID ring couple together strongly only under certain 
circumstances, specifically around particular values of 
the magnetostatic bias flux $> x . In this case, using the 
bias flux as a means to control the coupling, we have 
been able to reveal a whole range of interesting quantum 
phenomena. 

In previous workiiiiSii 3 . we dealt with the semi-classical 
problem of a monochromatic microwave field coupled to 
a SQUID ring containing a small capacitance weak link. 
In this paper we extend our theoretical description and 
treat both the ring and the field fully within a quantum 
mechanical framework. We demonstrate that the numer- 
ical results derived from this quantum model, in which 
the em field is initially in a coherent state, compare very 
well with those obtained using a semi-classical, Floquet 
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FIG. 1: A SQUID ring coupled inductively to a mode of an 
electromagnetic field. 



theory of a SQUID ring coupled to the field. In this 
there are obvious analogies to quantum optical interac- 
tions in few level atoms which apply to both pair con- 
densate and single electron systemsi^ii^iSiiLi&. In addi- 
tion, we note that SQUID rings have a strong sinusoidal 
non-linearity and it is the strength of this non-linearity, 
together with its periodic nature, that leads to the quite 
novel phenomena studied in this paper. This should be 
compared and contrasted with the large body of work on 
non-linear quantum systems in the context of quantum 
optics 19-20 ' 21 ' 22 - 23-24 where the non-linearity is usually a 
weak polynomial non-linearity. 



II. A SQUID RING COUPLED TO 
NON-CLASSICAL EM FIELD 



The Hamiltonian H t for our coupled system can be 
written as a sum of the energies for the field and the 
ring, together with an additional term for the interaction 
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energy, i.e. 



H t — H e + H s — H) 



Int- 



(1) 



where H e and H s are, respectively, the Hamiltonians for 
the field and the ring and -f/int is the interaction energy. 

We can write the Hamiltonian for the SQUID ring 
(weak link capacitance C s and ring inductance A s ), in 
the usual formii 



H„ 



2C S 



2A S 



hv cos 2tt 



$0 



(2) 



where $ s , the magnetic flux threading the ring, and Q s , 
the total charge across the weak link, are the conjugate 
variables for the system (with the imposed commuta- 
tion relation [<& S ,(5 S ] = ih), & x is the static (or quasi- 
static) external flux applied to the ring, hv/2 is the ma- 
trix element for pair tunnelling through the weak link 
(critical current I c — 2ev) and <f>o = h/2e. We note 
that with a characteristic frequency uj s = (l/VCsA s ) 
for the SQUID ring, there is a renormalized frequency 

1 related to the <I>^ term 



n. = 



in a Taylor expansion of the cosine in P|. Throughout 
the paper we use C s = 1 x 10~ 16 F, A s = 3 x 10~ 10 H 
and hv — 0.07$q/A s as typical circuit parameters for a 
SQUID ring in the quantum regime. 

The em field can be modelled in terms of a cavity mode 
using an equivalent circuit comprising a capacitance C e 
in parallel with an inductance A e , with a (parallel) re- 
sistance on resonance to define its quality factor. If we 
assume this resistance to be infinite we obtain a Hamil- 
tonian for the field in terms of the equivalent circuit flux 
and charge operators 



2C, 
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2A R 



(3) 



where $ e and Q e are, respectively, the magnetic flux and 
electric charge associated with the cavity. The field fre- 
quency is uj e = l/VC e A e . For the purposes of simplicity 
we use C e = C s throughout this paper and specify the 
frequency u> e in each example. We denote as |n) the 
eigenstates of H e . In our numerical work we use a trun- 
cated basis with n = 0, ...,7V, where N is taken to be 
much greater than the average number of photons in the 
system. 

The em cavity mode and the SQUID ring are coupled 
together inductively with a coupling energy given by 



#lnt = Y ($s 



(4) 



where fi is a coupling parameter linking the em field to 
the SQUID ring. 

We note that by introducing a unitary translation op- 
erator T = exp (—i§xQs/K) we can write the Hamilto- 
nian for the ring as 



O 2 $ 2 
Ht = T^H S T = + — ^ - fii/cos 



2C S 



2A.< 



2tt 



<S>s + <S>x 

<J>n 



(5) 
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FIG. 2: Energy eigenvalues versus tp x — & x /&o for an iso- 
lated SQUID ring. 



We also note that, invoking this unitary transformation, 
the interaction energy becomes H{ nt = -#-$ s $ e whilst 
the em field Hamiltonian remains unaffected. We denote 
as \a) the (flux-dependent) eigenstates of H' s . Again 
in our numerical work we use a truncated basis with 
a — 0, E, where E is taken to be much greater than 
the average energy level in which the SQUID operates. 

The first few eigenvalues (a = 0, 4) of H' s as functions 

of $> x /&o (— ip x ) are shown in figure 2. As can be seen, 
although all the eigenvalues are <E>o -periodic in & x , each 
displays a distinctive functional form in It will be- 
come apparent in the following discussion that these func- 
tional forms take on great importance in determining the 
behaviour of the coupled system at particular points in 
external bias flux. 

In describing the coupled system, we now intro- 
duce the dimensionless operators x e 
p e = ^l/C e huj e Q e , 



= ^C s uj s /h^ s and p s = 
^l/C s huj s Q s , together with the lowering and raising op- 



erators a s — -j= (x s + ip s ) 
ring and a e = ^= (x e + ip e ) 



t = J_ 



l 

t - 



x s — ip s ) for the 

the field. In terms of these operators the Hamiltonian 
H[ = T'HfT for the coupled system (see Q) can be 
rewritten in the form 



Hi = huj P [ aia e + - ) + huj s { aids 



1 



, 2n h 
hvcos | — — \j — x s + 2iripx 



$0 V CsLOs 
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As an illustrative example we show in figures 3 the 
computed, <p x -dependent eigenvalues of H[ for (a) u> e — 
uj s (with truncations N = E = 5) and (b) uj e = jqU s 
(with truncations N — 50 and E = 5 ). In these fig- 
ures the scaling is too small to reveal the lifting of the 
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FIG. 3: (a) Energy eigenvalues versus (p x of the SQUID ring 
Hamiltonian H s (thick lines) and the ring-field total Hamil- 
tonian H t (thin lines) with lj £ — uj s - The coupling constant 
(i = f/100. The inset shows an example (arrowed) of the lift- 
ing of the degeneracy of the ring-field levels when ^ / 0. (b) 
as in figure(a) but with w e = Jq^s- 



degeneracy at the crossing points by the nonzero cou- 
pling term (/i ^ 0). Again to illustrate, we show in the 
inset of figure 3(a), but at much higher resolution, one 
such computed crossing point. Here the splitting of the 
crossing energies is quite apparent, these being intimately 
connected to the functional form of the original SQUID 
ring eigenenergies. That such crossing points have been 
reported in experimental studies of Josephson weak link 
circuits, particularly SQUID rings, with concomitant su- 
perpositions of macroscopic states (Schrodinger cats), is 
further evidence for the underlying quantum mechanical 
nature of these systems§*2*2*iS. It is therefore timely to 
develop a full quantum treatment of SQUID ring-em field 
systems, which is the purpose of the paper. 

In the above we have studied the eigenproblcm 
H't = H n using a truncated basis. We use now 
these results to compute the evolution operator as 



CT(t)=5>»>exp(-^)<£ B 



(7) 
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FIG. 4: The time averaged energy levels versus ip x for (a) 
the ring ({H 3 ))and (b) the field ({H e )). The coupling constant 
(i = f/100 and u> e = ui s . At t — the electromagnetic field is 
assumed to be in the coherent state [a = iy/2) and the ring 
in the energy eigenstates: a = (blue), a = 1 (brown), a — 2 
(green) and a = 3 (red). 



density matrix p(0), we have calculated the density ma- 
trix pit) = V (t) p(0)W (t) at a later time t and the re- 
duced density matrices p e = Tr s (p) , p s = Tr e (p) . As 
a measure of the accuracy of the truncation approxima- 
tion we have also calculated the traces of all the density 
matrices that we use. In the limit of infinite order den- 
sity matrices the trace is equal to 1, while for truncated 
density matrices it should be very close to 1. In all our 
results the trace was greater than 0.99. Another test we 
performed was to increase the cutoff point from which we 
were able to ascertain that our truncation had negligible 
effect. 

We define the time averaged energy expectation values 
H' s and H e as (i ~ s, e) 



{(Hi)) 



1 f T 

lim - / Ti[ Pi (t)Hi] dt 



(8) 



Assuming that the system at t — is described by the 



where, computationally, we integrate from up to 
20, 000/o; s which we have found to be sufficient to ensure 
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FIG. 5: The time averaged energy levels versus tp x for (a) 
the ring ((_H' s })and (b) the field ({H e )). The coupling constant 
fi = I/I00 and lQuo e = uj s . At t — the electromagnetic field 
is assumed to be in the coherent state \a = ilOy/2) and the 
ring in the energy eigenstates: a = (solid), a = 1 (dashed). 



the convergence of the integral ijHJ for all the results pre- 
sented in this paper. In figure 4 we display the computed, 
time averaged, energy expectation values (normalized in 
units of fkv s ) of H' s [figure 4(a)] and H ' e [figure 4(b)]. 
These have been calculated over the range < ip x < 1 
for various values of a (= 0, 1, 2, 3), with li = 1/100 and 
uj e = cu s . In computing these results we have set the t = 
state as \a — iy2) e ® \a) 8 , where \a) e is a coherent state 
of the em field {a e \a) e = a\a) e ). As is apparent in fig- 
ures 4(a) and (b), for specific values of external bias flux 
ip x [namely those corresponding to the crossing points 
shown in figure 3(a)] there is a strong interaction between 
the em field and the SQUID ring. As is also apparent, 
this leads to an energy exchange between the components 
of the system. To demonstrate how the time averaged en- 
ergy levels for the ring and field depend on the ratio of ui a 
to oj e , we show in figure 5(a) and (b), respectively, these 
levels computed for a = Q, 1 and, again, \i = f/100 but 
now uj e = tOg/lO. In order for the energy of our initial 
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FIG. 6: The time averaged energy levels versus <p x for the ring 
({H s }) and the field ((H e )). The coupling constant /i = f/10 
and uj e — lo 3 . At t = the electromagnetic field is assumed to 
be in the coherent state [a = iy2) and the ring in the lowest 
eigenstate a — 0. 



state to be equal to that used in the previous example, 
here this state is chosen to be \a — zl0v / 2) e ® |c) s . As 
is to be expected, starting with a = 0, 1 eigenstates, the 
separation in ip x between the regions of strong coupling 
(energy exchange) are significantly reduced compared to 
those seen in figure 4. As a further example, we show 
in figure 6 the computed results for our coupled system 
taking, as in figure 4, uj e = lu s but now with stronger 
coupling (/i = 1/10). To make our results strictly com- 
parable with those of figure 4, we use the initial state 
\a = «v2) e ® W — 0)s- D ue to t ne stronger coupling we 
can see more regions in external bias flux where energy is 
exchanged between the two components of the system. In 
all three sets of results (figures 4, 5 and 6) there are peaks 
(both upwards and downwards) generated in the time av- 
eraged energies about specific values of <p x . These peak 
regions, where energy is exchanged between the field and 
the ring, correspond to quantum transitions in the ring 
and in all cases demonstrate strong coupling between the 
two oscillators in the system. 

To illustrate the quantum mechanical effects associ- 
ated with this coupling we take as an example the case 
of /i = 1/10, ui e = lo s (as in figure 6) and set <p x = 0.426 
(arrowed in figure 6) at which flux bias the coupling (and 
the energy exchange) between the ring and field is strong. 
We assume that at t = the em field is in the number 
state |1) (a\a e \l) e = l|l) e ) an d the ring is in the energy 
eigenstate a = (we stress that a are eigenstates of H 8 
and not of n s = a|a s ). In figure 7, with these values 
of u) s , \x and tp x , we show the computed expectation val- 



ues of the photon number 



Tr(/9 e aja e ) in the field, 



and (n s ) = Tr(p s a|a s ) in the ring, as functions of time. 
These results demonstrate that a strong exchange in en- 
ergy takes place quasi-periodically in time between the 
ring and the field, i.e. when the photon number expec- 
tation value in the SQUID ring increases that in the em 
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FIG. 7: The expectation values (n e ) and (n 3 ) as functions 
of time for tp x = 0.426, p = 1/10 and u) e = u> 3 . At t = the 
electromagnetic field is assumed to be in the number state 1 1} 
and the ring is in the energy eigenstate a — 0. 
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FIG. 8: The entropies as a function of time for the same 
system as in figure 7. 



field decreases, and vice versa. We note that in order 
to compare these predictions with experiment we would 
need to to measure the actual power level of the em field. 
We also note that with the ring-field coupling constant 
known this would allow us to estimate the em power im- 
pinging on the SQUID ring. 

As the system evolves in time its two components- 
oscillator mode and SQUID ring - become entangled 
quantum mechanically. In order to quantify this en- 
tanglement we use entropic quantities. For a two mode 
(field-ring) system e — s this entanglement can be quan- 
tified according to the expression25i2£i2LSL 

I es = S ( Pe ) + S ( Ps ) - S (p) (9) 

where S (p) is the von-Neumann entropy given by 

S(p) = -Tr[pln(p)} (10) 



the time dependent computed entropies S(p e ), S(p s ), 
and the entanglement entropy /, for the same system 
as in figure 7. From these results it is quite apparent 
that although at t = both the field and the ring are 
in a pure state, they both evolve into mixed states. Of 
course, since the time evolution is unitary the joint field- 
ring system is always in a pure state (S(p) = 0). The 
results presented in figure 8 do demonstrate that the sys- 
tem does become highly entangled over time although, 
as can be seen, at certain times it can disentangle again. 
There is no doubt that for the development of truly quan- 
tum technologies, for example, quantum computing and 
quantum communications, such states of entanglement 
are of great importance. We note in particular, that in 
some of these schemes the ability of the system to con- 
trol the entanglement (as in our case) is highly desirable. 
In principle, experimental verification of the entangle- 
ment between the two modes could be achieved through 
Bell type of inequalities. However, in the context of the 
present work their exact form will require further inves- 
tigation. 



III. A SQUID RING COUPLED TO A 
CLASSICAL EM FIELD 

In previous work we treated the em field classicalljiii 
and used the Hamiltonian 



H„ = 



01 

2C S 



($ s - [<S> X + tp e sm(u e t)]) 
2A, 



$ 



(11) 

and solved the corresponding time-dependent 
Schrodinger equation. Here tp e is taken to be the 
magnetic flux amplitude of the classical em field. It is 
of interest to compare the quasi-classical results derived 
via l|ll|) with the fully quantum results found above 
where the initial state of the em field is a coherent state. 
Due to the quasi-classical nature of the coherent state 
we expect some agreement between the fully quantum 
results and the quasi-classical results. To furnish an 
example to compare with these quantum results, we have 
computed the time averaged ring energy expectation 
values for the Floquet states (eigenvalues of the evolution 
operator after one period of microwave evolution) as a 
function of (p x using the same value of microwave field 
amplitude (<p e = 0.41/i<f>o) as in figure 4(a). Our results 
are presented in figure 9. Within the computational 
accuracy available, and given that we are dealing with 
two different regimes of the coupled system, it is clear 
that the principal transition region features match in 
both models, even though the amplitudes may not be 
the same. 



IV. DISCUSSION 



and the entanglement entropy is positive or zero (sub- 
additivity property of the entropy). In figure 8 we show 



We have studied the coupling of a SQUID ring to a 
single mode em field at frequencies in the sub-THz to 
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FIG. 9: Time-dependent Schrodinger equation (Eq 12) cal- 
culation of the time averaged ring energy levels against <& x 
with n — 1/100 , uj e — w s and tp e — 0.41/^$o Here the electro- 
magnetic field is treated classically and the time-dependent 
Schrodinger equation (12) is solved. The results of this figure 
should be compared and contrasted with those of figure 4a. 



THz range where, in general, ksT < huj s ,fiw e so that 
the system behaves quantum mechanically. In this we 
have been strongly influenced by recent discussions in 
the literature on routes to quantum computing using 
Josephson devicesi*2ii4 and on even more recent pub- 
lication of experimental data on superposition states in 
SQUID rings&iil. In the paper we have shown that our 
results [for example, figure 4(a)] compare well with previ- 
ous semi-classical (Floquet method) computations made 
by us (figure 9 and ref 11) and have expanded this work to 
calculate explicitly the photon number and entanglement 
states in the coupled ring-field system. We note the rela- 
tion between the results presented in this paper and pre- 
vious work on em environments in thermal equilibrium 
with Josephson (e.g. SQUID) circuits 29 i 30 i 31 i 32 i 33 i 34 i 35 . 
In our work we have assumed that the external em field 
is in a particular quantum state and our results depend 
on this state. As an example we have considered the 
em field at t — to be in a coherent state. However, 
the calculations can easily be repeated for another initial 
state of the field. Each initial state will, of course, yield 
different results. 

In addition to the em field a magnetostatic flux also 
threads the SQUID ring. In this paper we have demon- 
strated that for certain values of this flux strong coupling 
develops at which point (s) large amounts of energy are 
exchanged between the ring and the field. Future ex- 
perimental probing of these energy exchanges, which is 
considered again in the discussion section, would clearly 
be of great interest. We have also demonstrated that en- 
tanglement between the ring and field modes arises as 
a natural consequence of the full quantum mechanics of 
the system. As we have pointed out (section II, above), 
experimental verification of such entanglements will re- 
quire further work on the Bell inequalities related to these 



entanglements. 

In our calculations we have neglected dissipation and 
have calculated the time evolution of the system us- 
ing the equation d t p — — (i/h) [H, p]. A more realistic 
calculation to take into account dissipation due to the 
external environment^ can be made with the equation 
d t p = — (i/h) [H, p] + f(p) where the f(p) are dissipative 
terms. Numerical work to include these terms is currently 
in progress. 

Of general interest to experimentalists working on 
(time dependent) superposition states in SQUID ring de- 
vices is the problem of determining the actual em power 
(or number of photons in each state of the em field) cou- 
pled to the ring. Together with the frequency of the 
em field, and the original eigenenergies of the ring, this 
is required in order to compute the ring crossing point 
splittings. In principle, this problem can be overcome 
through the kind of analysis we have undertaken in this 
paper, whether it be for classical em fieldsAi or photon 
states interacting with a quantum mechanical SQUID 
ring. As we have shown, it is possible to determine these 
power levels accurately through following the reactive fre- 
quency shift of a SQUID ring-classical (radio frequency) 
resonator system, driven by an external em field, when 
the ring remains adiabatically in its quantum mechanical 
ground state. However, where em frequencies and/or am- 
plitudes are high enough (as in this paper), so that (time 
dependent) superpositions of low lying energy eigenstates 
of the ring are generated, the problem becomes very 
much more difficult theoretically. Nevertheless, there ap- 
pear to be several routes to resolving these difficulties, 
as indicated by some of our recent investigations of non- 
adiabatic processes in em-driven SQUID rings We 
note that at sufficiently high em frequencies/amplitudes 
multiphoton absorption and emission processes will oc- 
cur between the components of the coupled system. This 
may complicate the interpretation of experimental data 
and will be a topic of further theoretical investigation by 
us. For the future, we also note that it may be possible to 
extend these experimental and theoretical techniques to 
probe the details of energy exchange and entanglement 
of the system presented in this paper. 

There now exists a clearly defined need to create THz 
technology 3 ^ for a range of applications including modern 
communications. To date this technology, based on quan- 
tum processes on the small scale, functions classically at 
the device level. A primary purpose of this paper has 
been to demonstrate that at THz frequencies, and rea- 
sonable operating temperatures (~4if), this technology 
could be made fully quantum mechanical in nature, i.e. 
at high enough SQUID ring and em oscillator frequencies 
both can be treated as macroscopic quantum objects, 
irrespective of any deeper description of the supercon- 
ducting condensate in SQUID rings 40 - 41 . This would 
point to a great richness of potential applications. For 
example, in the context of the results presented here, our 
investigations may prove useful for the development of 
frequency converters up to THz frequencies and beyond. 



7 



More generally these results, and the theoretical descrip- itude to Professors C.H. van de Wal and A.Sobolev for 
tions underlying them, may find use in the emerging fields interesting and informative discussions, 
of quantum technology and quantum computatio n 1 ! 2 ^ 4 . 
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